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Abstract 

In this work we study the homogenization for infinitesimal disloca- 
tion based gradient viscoplasticity with linear kinematic hardening and 
general non-associative monotone plastic flows. The constitutive equa- 
tions in the models we study are assumed to be only of monotone type. 
Based on the generalized version of Korn's inequality for incompatible 
tensor fields (the non-symmetric plastic distortion) due to Neff et al., we 
derive uniform estimates for the solutions of quasistatic initial-boundary 
value problems under consideration and then using an unfolding operator 
technique and a monotone operator method we obtain the homogenized 
system of equations. 
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1 Introduction 

We study the homogenization of quasistatic initial-boundary value problems 
arising in gradient viscoplasticity. The models we study use rate-dependent 
constitutive equations with internal variables to describe the deformation be- 
haviour of metals at infinitesimally small strain. 

Our focus is on a phenomenological model on the macroscale not including 
the case of single crystal plasticity. From a mathematical point of view, the maze 
of equations, slip systems and physical mechanisms in single crystal plasticity 
is only obscuring the mathematical structure of the problem. 
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Our model has been first presented in [37]. It is inspired by the early work 
of Menzel and Steinmann |33| . Contrary to more classical strain gradient ap- 
proaches, the model features from the outset a non-symmetric plastic distortion 
field p G M.^ 0, a dislocation based energy storage based solely on | Curlpj (and 
not Vp) and therefore second gradients of the plastic distortion in the form of 
CurlCurlp acting as dislocation based kinematical backstresses. We only con- 
sider energetic length scale effects and not higher gradients in the dissipation. 

Uniqueness of classical solutions in the subdifferential case (associated plas- 
ticity) for rate-independent and rate-dependent formulations is shown in |36) . 
The existence question for the rate-independent model in terms of a weak refor- 
mulation is addressed in [37] . The rate-independent model with isotropic hard- 
ening is treated in |19] . First numerical results for a simplified rate-independent 
irrotational formulation (no plastic spin, symmetric plastic distortion p) are pre- 
sented in [H] . In [521 EO] well-posedness for a rate-independent model of Gurtin 
and Anand [24] is shown under the decisive assumption that the plastic distor- 
tion is symmetric (the irrotational case), in which case one may really speak of 
a strain gradient plasticity model, since the full gradient acts on the symmetric 
plastic strain. 

In order to appreciate the simplicity and elegance of our model we sketch 
some of its ingredients. First, as is usual in infinitesimal plasticity theory, 
we split the total displacement gradient into non-symmetric elastic and non- 
symmetric plastic distortions 

Vu ^ e -\- p . 

For invariance reasons, the elastic energy contribution may only depend on the 
elastic strains syme = sym(Vu — p). While p is non-symmetric, a distinguishing 
feature of our model is that, similar to classical approaches, only the symmetric 
part Ep := symp of the plastic distortion appears in the local Cauchy stress ct, 
while only the higher order stresses are non-symmetric. The reason for this is 
that we assume that p has to obey the same transformation behavior as Vu 
does, and thus the energy storage due to kinematical hardening should depend 
only on the plastic strains symp. For more on the basic invariance questions 
related to this issue dictating this type of behaviour, see [54j[35]. We assume as 
well plastic incompressibility trp = 0, as is usual. 

The thermodynamic potential of our model can therefore be written as 

/ (c[a;](sym(Vu-p))(sym(Vu-p)) (1) 

clastic energy 

+ ^^|devsymp|2-f ^jCurlpp + 

^ V ^ V ^ external volume forces 

kinematieal hardening disloeation storage 

The positive definite elasticity tensor C is able to represent the elastic anisotropy 
of the material. The evolution equations for the plastic distortion p are taken 
such that the stored energy is non-increasing along trajectories of p at frozen 
displacement u, see [37]. This ensures the validity of the second law of thermo- 
dynamics in the form of the reduced dissipation inequality. 

For the reduced dissipation inequality we fix u in time and consider the time 
derivative of the free energy (and taking into account that Curl is a self-adjoint 
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operator provided that the appropriate boundary conditions are specified), we 
have 

d f 

- W{Wu{ta) ~ Pit), Pit), Curlp(i)) dx 

= / DiW ■ i-dtp)+ D2W -dtp + D3W -Cml dtp dx 
Jn 

= -( iDiW ~D2W -Cni\D:iW)-dtpdx. 
Jn 

Choosing dtp E giDiW — D2W — Curli^aW^) with a monotone function g we 
obtain the reduced dissipation inequality 

- / Wi\/uito)-pit),pit),Cm\pit))dx<0. 
J o 

Adapted to our situation, the plastic flow has the form 

dtp e gi<T — Ci [x] dev symp — C2 Curl Curlp) , (2) 

where a — C[a;]sym(Vu — p) is the elastic symmetric Cauchy stress of the 
material and g is a multivalued monotone flow function which is not necessary 
the subdifferential of a convex plastic potential (associative plasticity). In this 
generality, our formulation comprises certain non-associative plastic flows in 
which the yield condition and the flow direction are independent and governed 
by distinct functions. Moreover, the flow function g is supposed to induce a 
rate-dependent response as all materials are, in reality, rate-dependent. 

Clearly, in the absence of energetic length scale effects (C2 — 0), the Curl Curlp- 
term is absent. In general we assume that g maps symmetric tensors to symmet- 
ric tensors. Thus, for C2 = the plastic distortion remains always symmetric 
and the model reduces to a classical plasticity model. Therefore, the energetic 
length scale is solely responsible for the plastic spin in the model. The appear- 
ance of the Curl Curl p-term in the argument of g is clear: the argument of g 
consists of the Eshelby-stress tensor E driving the plastic evolution, see [37] . 

Regarding the boundary conditions necessary for the formulation of the 
higher order theory we assume that the boundary is a perfect conductor, this 
means that the tangential component of p vanishes on d^l. In the context of 
dislocation dynamics these conditions express the requirement that there is no 
flux of the Burgers vector across a hard boundary. Gurtin [25 introduces the 
following different types of boundary conditions for the plastic distortior0 

dtp X r^lrhard — "micro-hard" (perfect conductor) 

9tp|rha.d=0 "hard-slip" (3) 
Curlp X n|ri,^^d = "micro- free" . 

We specify a sufficient condition for the micro-hard boundary condition, namely 

and assume for simplicity Fhard = dQ. This is the correct boundary condition 
for tensor fields in ff(Curl) which admits tangential traces. 



-"^Here, v X n with v G and n G R'' denotes a row by column operation. 
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We combine this with a new inequahty extending Korn's inequahty to in- 
compatible tensor fields, namely 

3C = C{n) > Vp e i?(Curl) ; p x n|r,„<j - : (4) 
lbllL2(n) < C(fi) II symj3||i2(Q) + [| Curlp||L2(o) ) . 

plastic distortion plastic strain dislocation density 

Here, Fhard C with full two-dimensional surface measure and the domain 
VI needs to be sliceable, i.e. cuttable into finitely many simply connected 
subdomains with Lipschitz boundaries. This inequality has been derived in |38l 
[551 SO] and is precisely motivated by the well-posedness question for our model 
[57) . The inequality (|H) expresses the fact that controlling the plastic strain 
symp and the dislocation density Curlp in L^(il) gives a control of the plastic 
distortion p in L^(f2) provided the correct boundary conditions are specified: 
namely the micro-hard boundary condition. Since we assume that tr(p) = 
(plastic incompressibility) the quadratic terms in the thermodynamic potential 
provide a control of the right hand side in Q. 

It is worthy to note that with g only monotone and not necessarily a subdif- 
ferential the powerful energetic solution concept [321 I5T| cannot be applied. 
In our model we face the combined challenge of a gradient plasticity model 
based on the dislocation density tensor Curlp involving the plastic spin, a gen- 
eral non-associative monotone flow-rule and a rate-dependent response. 



Setting of the homogenization problem. Let C M"^ be an open bounded 
set, the set of material points of the solid body, with a C^-boundary and y C M"^ 
be a set having the paving property with respect to a basis [hi , 62 , 63) defining the 
periods, a reference cell. By Tg we denote a positive number (time of existence), 
which can be chosen arbitrarily large, and for < < < Te 

nt^Vlx (0,i). 

The sets, and denote the sets of all 3 x 3-matrices and of all symmetric 
3 X 3-matrices, respectively. Let s[(3) be the set of all traceless 3 x 3-matrices, 
i.e. 

s[(3) = {w e I tri) = 0}. 

Unknown in our small strain formulation are the displacement Ujj{x,t) G 
of the material point x at time t and the non-symmetric infinitesimal plastic 
distortion Pri{x, t) £ s[(3). 

The model equations of the problem are 

- div^o-^(a;,i) = b{x,t), (5) 
<Jnix,t) = C[x/r]]{sym(yxUrjix,t) - pr,{x,t))), (6) 

dtp,{x,t) e g{x/rj,j:^^{x,t)), Ef = E:,';; Elt, + ^curi,,' (7) 
= a„ El"^ = -Ci [x/rj] dev symp^, j:'Zi,r, = "^2 Curl Curlp^ , 



which must be satisfied in x [0,Te). Here, C2 > is a given material con- 
stant independent of 77 and Ejj" is the infinitesimal Eshelby stress tensor driving 
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the evolution of the plastic distortion and 77 is a scaling parameter of the 
microstructure. The initial condition and Dirichlet boundary condition are 

p^{x,o) = xen, (8) 

Pr^{x,t) X n{x) = 0, {x,t) e dnx [0,Te), (9) 

u^{x,t) = 0, {x,t) ednx[0,Te), (10) 

where n is a normal vector on the boundary For simplicity we consider 

only homogeneous boundary condition and we assume that the cell of periodicity 
is given by y = [0, 1)^. Then, we assume that Ci : F — M, a given material 
function, is measurable, periodic with the periodicity cell Y and satisfies the 
inequality 

Ci[y]>ai >0 

for all y G y and some positive constant ai. For every y £ Y the elasticity 
tensor C[?;] : 5'^ 5^ is a linear, symmetric, positive definite mapping and the 
mapping y i— > C[j/] : — > 5^ is measurable and periodic with the same periodic- 
ity cell Y. Due to the above assumption (Ci > 0), the classical linear kinematic 
hardening is included in the model. Here, the nonlocal backstress contribution 
is given by the dislocation density motivated term ^ = —C2 Curl Curl 
together with corresponding Neumann conditions. 

For the model we require that the nonlinear constitutive mapping v t-^ 
g{y,v) : — >■ 2^'('^^ is monotone for all y &Y, i.e. it satisfies 

< iv,-V2)-iv*,-v;), (11) 

for all Vi € M^, v* G g{y, Vi), i — 1,2 and all y £ Y. We also require that 

Oe5(2/,0), a.e.yeY. (12) 

The mapping y ^ g{y, •) : 2^^'-^^ is periodic with the same periodicity cell 

Y. Given are the volume force b{x, t) e M.^ and the initial datum p^'^^x) G s[(3). 

Remark 1.1. It is well known that classical viscoplasticity (without gradient ef- 
fects) gives rise to a well-posed problem. We extend this result to our formulation 
of rate- dependent gradient plasticity. The presence of the classical linear kine- 
matic hardening in our model is related to Ci > whereas the presence of the 
nonlocal gradient term is always related to C2 > 0. 

The development of the homogenization theory for the quasi-static initial 
boundary value problem of monotone type in the classical elasto/visco-plasticity 
introduced by Alber in [2] has started with the work [31 , where the homogenized 
system of equations has been derived using the formal asymptotic ansatz. In 
the following work [3] Alber justified the formal asymptotic ansatz for the case 
of positive definite free energjlf], employing the energy method of Murat-Tartar, 
yet only for local smooth solutions of the homogenized problem. It is shown 
there that the solutions of elasto/visco-plasticity problems can be approximated 
in the L^(f2)— norm by the smooth functions constructed from the solutions of 
the homogenized problem. Later in [42 , under the assumption that the free 
energy is positive definite, it is proved that the difference of the solutions of 

^Here, v X n with v € and n € R'^ denotes a row by column operation. 

^Positive definite energy corresponds to linear kinematic hardening behavior of materials. 
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the microscopic problem and the solutions constructed from the homogenized 
problem, which both need not be smooth, tends to zero in the L^(ri x F)— norm, 
where Y is the periodicity cell. Based on the results obtained in [32], in [S] the 
convergence in L^(0 x Y) is replaced by convergence in L^{U). In the meantime, 
for the rate-independent problems in plasticity similar results are obtained in 
[34) using the unfolding operator method (see Section [3]) and methods of ener- 
getic solutions due to Mielke. For special rate-dependent models of monotone 
type, namely for rate-dependent generalized standard materials, the two-scale 
convergence of the solutions of the microscopic problem to the solutions of the 
homogenized problem has been shown in [56l [57] . The homogenization of the 
Prandtl-Reuss model is performed in [52j [57]. In [43] the author considered 
the rate-dependent problems of monotone type with constitutive functions g, 
which need not be subdifferentials, but which belong to the class of functions 
M(ri, M^,q, a, m) introduced in Section[5| Using the unfolding operator method 
and in particular the homogenization methods developed in [T^ , for this class of 
functions the homogenized equations for the viscoplactic problems of monotone 
type are obtained in [15] . 

In the present work the construction of the homogenization theory for the 
initial boundary value problem ([S]) - pH]) is based on the existence result derived 
in [3S] (see Theorem I5.6p and on the homogenization techniques developed in 
[43] for classical viscoplasticity of monotone type. The existence result in [45] 
extends the well-posedness for infinitesimal dislocation based gradient viscoplas- 
ticity with linear kinematic hardening from the subdifferential case (see [H]) to 
general non-associative monotone plastic flows for sliceable domains. In this 
work we also assume that the domain f2 is sliceable and that the monotone 
function g : M? x M'^ 2^'^^) belongs to the class M(r2, A^^, g, a, m). For 
sliceable domains £7, based on the inequality ([4]), we are able to derive then uni- 
form estimates for the solutions of ([S]) - (flUl) in Lemma 15.81 Using the uniform 
estimates for the solutions of ([5]) - (ITU|) . the unfolding operator method and 
the homogenization techniques developed in [THI |33] , for the class of functions 
M[(ri, q, a, m) we obtain easily the homogenized equations for the original 
problem under consideration (see Theorem 15. 7|) . To the best our knowledge 
this is the first homogenization result obtained for the problem ([SJ - pUj) . We 
note that similar homogenization results for the strain-gradient model of Fleck 
and Willis [10] are derived in [^T] [231 [2Z] using the unfolding method together 
with the F-convergence method in the rate-independent setting. In [3T] the 
authors, based on the assumption that the model under consideration is of rate- 
independent type, are able to treat the case when C2 is a F-periodic function 
as well. In the rate-independent setting this is possible due to the fact that the 
whole system ([SJ - (jlOp can be rewritten as a standard variational inequality 
(see [5S] ) and then the subsequant usage of the techniques of the convex anal- 
ysis enable the passage to the limit in the model equations. Contrary to this, 
in the rate-independent case this reduction to a single variational inequality is 
not possible and one is forced to use the monotonicity argument to study the 
asymptotic behavior of the third term Sj,'"^! ,j in ([7]). 

Notation. Suppose that f2 is a bounded domain with C°°-boundary dO.. 
Throughout the whole work we choose the numbers g, q* satisfying the following 
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conditions 

l<q,q*<oo and l/q+l/q* = 1, 

and I • I denotes a norm in M*^. Moreover, the following notations are used in 
this work. The space Ty"''^(r2, M*^) with q G [l,oo] consists of all functions in 
L''( fl, K''") with weak derivatives in L'^{fl, M'^) up to order m. If m is not integer, 
then M''') denotes the corresponding Sobolev-Slobodecki space. We set 

if'"(fi,R'=) = W"'^^{n,R''). The norm in is denoted by || • |U,5,n 

(II • 11? II ■ Wo.qS-i)- The operator Tq defined by 

denotes the usual trace operator. The space W^''^{ft,B}^) with q G [l,oo] con- 
sists of all functions v in VF™'*(f2, R'^) with Tqv = 0. One can define the bilinear 
form on the product space L'^{il, M^)xL'^' {il,, M^) by 

(C,C)a= / a^)-C{x)dx. 
Jn 

The space 

^Curi(^'-^^) = e L'J{n,M^) I Curlw G L'J{n,M^)} 
is a Banach space with respect to the norm 

lk|l9,Curl= |k||, + ||Curll>||,. 

By _ff(Curl) we denote the space of measurable functions in L'^^^i{i}, M^), i.e. 
i?(Curl) = L^j^j.[(r2, The well known result on the generalized trace oper- 

ator can be easily adopted to the functions with values in (see [SSI Section 
II. 1.2]). Then, according to this result, there is Si bouiid.Gd operator r^^, on 

with 

r„w = v X n\g^ if V G C^{Q,,M^), 
where X* denotes the dual of a Banach space X. Next, 

L'curi.oi^^M') ^{we Ll^^,{n,M') I r„(«;) = 0}. 

Let us define spaces V'i{n,M^) and X'i{n,M^) by 

V''{n,M^) = {v G L'^{n,M^) I divw,Curlw G Li{n,M'^),rr,v = 0}, 

X'^{n,M^) = {ve L'^{n,M^) I divv,Curlv G L^in, M^),rQV = 0}, 
which are Banach spaces with respect to the norm 

ll«lk<!(lklUO = lkll« + II Curlwllg + II divvllg. 

According to [30l Theorem 20 the spaces V''{n,M^) and X'i{n,M^) are con- 
tinuously imbedded into W^^i{n,M^). We define V^{n,M^) and X^{n,M^) 

by 

V^{n,R^) {v G F9(1),R3) I divv = 0}, 

^This theorem has to be apphed to each row of a function with values in Ai^ to obtain the 
desired result. 
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x^{n,R^) := {v e x'i{n,R^) \ diw = o}, 

and denote by Vl^^^{V,, R^) and X^^^{i}, M.^) the L''-spaces of harmonic functions 
on n, as 

Vl^in^W") {v e V,\n,R^) I Curli, = 0}, 

Xl^^{n,W') := {v € X^{n,W') \ Cmlv = 0}, 

Then the spaces V/l^^{n,M.^) and X^^^(rj,R3) for every fixed q, I < q < oo, 
coincides with the spaces Vhar{^,^'^) and K"^) given by 

Vhari^,R^) = {v e C°°{n,R^) I divw = 0,Curlw = with w • n = on 90}, 

Xhari^,^^) ^{v e C°°(n,R3) I divw = 0,Curlw = with w X n = on dfl}, 

respectively (see [301 Theorem 2.1(1)]). The spaces Vhari^,^^) and X?,ar(0, M^) 
are finite dimensional vector spaces ( |30[ Theorem 1]). 
We also define the space Z^^^^{n, M^) by 

ZLA^^M') = e Ll^nifli^^M') I Curl Curl z; e L''{n,M'')}, 
which is a Banach space with respect to the norm 

II^II^SuH ll'^II^^Curi + II Curl Curl wllg. 

The space W^;?{Y, M'^) denotes the Banach space of y-periodic functions in 
W^f{R'',m'') equipped with the VF'"'9(y,M'=)-norni. 

For functions v defined on J7 x [0, oo) we denote by v{t) the mapping x i— > 
v{x, t), which is defined on fl. The space ^'(0, Te] X) denotes the Banach space 
of all Bochner-measurable functions u : [0,Te) — > X such that t i— >■ ||u(t)||^ is 
integrable on [0, Te). Finally, we frequently use the spaces VF'"'^(0, Tg; X), which 
consist of Bochner measurable functions having g-integrable weak derivatives up 
to order m. 

2 Maximal monotone operators 

In this section we recall some basics about monotone and maximal monotone 
operators. For more details see [8ll28l|48], for example. 

Let y be a reflexive Banach space with the norm || • ||, V* be its dual space 
with the norm || • ||*. The brackets (•, •} denotes the dual pairing between V and 
V*. Under V we shall always mean a reflexive Banach space throughout this 
section. For a multivalued mapping A:V ^2^' the sets 

D{A) = {veV \Av^9} 

and 

GrA ^ {[v,v*] eV xV* \ v e D{A), v* e Av] 
are called the effective domain and the graph of A, respectively. 

Definition 2.1. A mapping A : V ^ 2^ is called monotone if and only if the 
inequality holds 

{v* -u*,v-u)>0 V [v, u*], [m, u*] e GrA. 
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A monotone mapping A : V 2^ is called maximal monotone iff the 
inequality 

{v*-u\v-u)>0 y[u,u*]£GrA 

implies [v,v*] G GrA. 

A mapping A : V 2^ is called generalized pseudomonotone iff the set Av 
is closed, convex and bounded for all v G D{A) and for every pair of sequences 
{vn} and {v!^} such that w* £ Avn, Vn vq, u* Wq G V* and 

limsup (u*,w„ - Wo) < 0, 

n— >oo 

we have that [wo,Wq] S GrA and (-y*,w„) {vq,vo). 

A mapping A : V 2^ is called strongly coercive iff either D{A) is bounded 
or D{A) is unbounded and the condition 

77— r. > +00 as \\v\\ — > cx), \v,v e GrA, 

\M 

is satisfied for each w G D{A). 

It is well known (051 P- 105]) that if A is a maximal monotone operator, 
then for any v G D{A) the image Av is a closed convex subset of V* and the 
graph GrA is demi-closedH A maximal monotone operator is also generalized 
pseudomonotone (see [HlUHlllH])- 

Remark 2.2. We recall that the subdifFerential of a lower semi-continuous and 
convex function is maximal monotone (see [49, Theorem 2.25]). 

Definition 2.3. The duality mapping J : V 2^ is defined by 
J{v)^{v*eV* I {v*,v)^\\vr^\\v*\\l } 

for all v ^ V . 

Without loss of generality (due to Asplund's theorem) we can assume that 
both V and V* are strictly convex, i.e. that the unit ball in the corresponding 
space is strictly convex. In virtue of [8, Theorem II. 1.2], the equation 

J{vx -v) + XAv\ 3 

has a solution vx G D{A) for every v € V and A > if A is maximal monotone. 
The solution is unique (see jB] p. 41]). 

Definition 2.4. Setting 

vx — jxv and Axv = —X~^J{vx — v) 

we define two .single valued operators: the Yosida approximation Ax : V ^ V* 
and the resolvent fx - V ^ D{A) with D{Ax) ^ D{j^) ^ V . 

By the definition, one immediately sees that Axv G A(^j^v^. For the main 
properties of the Yosida approximation we refer to [H [251 HE] and mention 
only that both are continuous operators and that Ax is bounded and maximal 
monotone. 

set A S V X V is demi-closed if d„ converges strongly to vq in V and ii* converges 
weakly to Vg in V* (or d„ converges weakly to vq in V and f* converges strongly to Vg in 
V*) and [vn,v^] g GrA, then [v,v*] € GrA 
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Convergence of maximal monotone graphs In the presentation of the 
next subsections we follow the work \i6\, where the reader can also find the 
proofs of the results mentioned here. 

The derivation of the homogenized equations for the initial boundary value 
problem ([5]) - (jlOl) is based on the notion of the convergence of the graphs of 
maximal monotone operators. 

According to Brezis [TD] and Attouch [7] , the convergence of the graphs of 
maximal monotone operators is defined as follows. 

Definition 2.5. Let A"', A : V 2^ he maximal monotone operators. The 
sequence v4" converges to A as n oo, (A^^ A), ij for every G GrA 

there exists a sequence [w„,w*] G GrA^ such that [w„,w*] — >■ strongly in 

V xV* as n ^ oo. 

Obviously, if A^ and A are everywhere defined, continuous and monotone, 
then the pointwise convergence, i.e. if for every w S y, ^"(w) — ^ ^(1^), implies 
the convergence of the graphs. The converse is true in finite-dimensional spaces. 

The next theorem is the main mathematical tool in the derivation of the 
homogenized equations for the problem ([5]) - (TlUl) . 

Theorem 2.6. Let A^, A : V ~^ 2^ he maximal monotone operators, and 
let [w„,w*] G GrA" and [v,v*] e V x V* . If, as n 00, A" A, Vn vq, 
w* ^ G V* and 

limsup«,w„) < {vq,vo) , (13) 

n— f 00 

then [vo,Vq] G GrA and 

liminf = {vo,vo) . 

n— ^oo 

Proof See W, Theorem 2.8]. □ 

The convergence of the graphs of multi- valued maximal monotone operators 
can be equivalently stated in term of the pointwise convergence of the corre- 
sponding single- valued Yosida approximations and resolvents as the following 
result shows. 

Theorem 2.7. Let A"', A : V ^ 2^ be maximal monotone operators and 
A > 0. The following statements are equivalent: 

(a) A" A as n 00; 

(h) for every v V, w — > j^v as n 00; 

(c) for every v £V , A^ \v — ?> A\v as n ^ 00; 

(d) A^l ^ A\ as n 00. 

Moreover, the convergences v j^v and A^^xv ~> Axv are uniform on 
strongly compact suhsets of V . 

Proof See [M Theorem 2.9]. □ 
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Measurability of multi-valued mappings. In this subsection we present 
briefly some facts about measurable multi-valued mappings. We assume that V, 
and hence V*, is separable and denote the set of maximal monotone operators 
from V to V* by m{V xV*). Further, let (5, T,{S), n) be a a-finite /i-complete 
measurable space. 

Definition 2.8. A junction A : S ^ 9Jl(T^ x V*) is measurable iff for every 
open set U x V* (respectively closed set, Borel set, open hall, closed ball), 

{x(^S\A{x)^U^%} 

is measurable in S. 

The next result states that the notion of measurability for maximal monotone 
mappins can be equivalently defined in terms of the measurability for appropri- 
ate single- valued mappings. 

Proposition 2.9. Let A : S ^ M{V x V*), let \ > and let E be dense in V . 
The following are equivalent: 

(a) A is measurable; 

(b) for every v ^ E , x i-^ Jx^^'^'^ measurable; 

(c) V Cz E, a;i— > Ax{x)v is measurable. 

Proof See [M Proposition 2.11]. □ 
For further reading on measurable multi- valued mappings we refer the reader 

to [misHiiiz]- 

Canonical extensions of maximal monotone operators. Given a map- 
ping A : S -n- dJl{V X V*), one can define a monotone graph from LP{S, V) to 
L«(S',y*), where l/p+ l/g= 1, as follows: 

Definition 2.10. Let A : S dJl{V x V*), the canonical extension of A from 
LP{S, V) to L^iS, V*), where l/p+l/q = l, is defined by: 

GrA^{[v,v*] £ LP{S,V)xL'^{S,V*) \ [v{x),v*{x)] E GrA{x) for a.e. x G S}. 

Monotonicity of ^ defined in Definition 12 . 101 is obvious, while its maximality 
follows from the next proposition. 

Proposition 2.11. Let A : S ^ m{V x V*) be measurable. If GrA ^ 0, then 
A is maximal monotone. 

Proof See [M Proposition 2.13]. □ 

We have to point out here that the maximality of A{x) for almost every 
X E S does not imply the maximality of A as the latter can be empty (^|16)): 
S = (0, 1), and GrA{x) = {[v,v*] E R x R \ v* = x"!/-?}. 

For given mappings A, A" : S — J> D}l{V x V*) and their canonical extensions 
A, A"", one can ask whether the pointwise convergence ^ A{x) implies 

the convergence of the graphs of the corresponding canonical extensions >^ 
A. The answer is given by the next theorem. 
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Theorem 2.12. Let A,A'^ : S ^ M{V x V*) be measurable. Assume 

(a) for almost every x G S, j4"(x) >— > A{x) as n ^ oo, 

(b) A and are maximal monotone, 

(c) there exists [«„,/?„] G Gr^" and [a,p\ S Lp{S,V) x L9(S',T/*) such that 
[a,f3] [a, /3] strongly in LP{S, V) x L''{S,V*) as oo, 

then A"- ^ A. 

Proof. See [161 Proposition 2.16]. □ 

We note that assumption (c) in Theorem 12 . 1 2 1 can not be dropped in virtue 
of Remark 2.16 in [TB]. 



3 The periodic unfolding 

The derivation of the homogenized problem for ([5]) - (jlOp is based on the periodic 
unfolding operator method introduced by Cioranescu, Damlamian and Griso 
[12| . For the reader unfamiliar with this method we recall here some properties 
of this operator. The proofs of all results mentioned here as well as examples 
of applications of the method can be found in [T^l UHl [H] and in the literature 
cited there. 

Let SI C M'^ be an open set and Y = [0,1)'^. Let (61,62,63) denote the 
standard basis in R^. For 2; G M"^, [z]y denotes a linear combination '^j^j 
with {di,d2, da} G Z such that z — [z]y belongs to Y, and set 

{z}y ■■= z - [z]y e Y ugM^ 

Then, for each x G M"^, one has 





X 








+ y) 




.V. 


Y / 



We use the following notations: 

E„ = G I rK? + r) c n}, = mt | y (r,e + ryF) I , A^^n\ n^. 

The set Cl,j is the largest union of 77(^ + y) cells G Z'^) included in il, while 
Arj is the subset of ft containing the parts from ri{£^ +Y) cells intersecting the 
boundary dfl. 

Definition 3.1. Let Y be a reference cell, rj be a positive number and a map 
V : il ^ M.'' . The unfolding operator 7^(w) : x Y ^ R'^ is defined by 




a.e. {x, y) € rtrj X Y, 
a.e. {x, y) G X Y. 

The next results are straightforward from Definition 13.11 
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Proposition 3.2. For q G [l,oo[, the operator Tri is linear and continuous from 
L«(f7,R'=) to L'J{n X Y^M.^). For every <p in L^{n,R'') one has 

(^) W\ inyY^-ni4>){x,v)dxdy = /^^ (t>{x)dx, 

W M^xy \%{4'){x,y)\dxdy < ^ \(t)ix)\dx, 

M In^ (k{x)dx - 1^ /j^xy '^('?^)(^' y)dxdy < J^^ \(j}{x)\dx, 

(d) \\%{<t>)\\p,n.Y - \Y\''n^I0, < \Y\'^m\,. 

Proof. See [13^ Proposition 2.5]. □ 
Obviously, if 0^ G L^{Q,,R'') satisfies 

\cl)^{x)\dx ^ 0, (14) 



then 

1 



Y 



(l}ri{x)dx - — T,-,{(j}n){x,y)dxdy ^ Q. 



If a sequence satisfies (|14p . we sliall write 

T„ 1 /■ 

<l)r^{x)dx ~ — %{(t)r^){x,y)dxdy. 
\^ I JfJxy 

Proposition 3.3. Lef g belong to [l,oo[. 

(a) For any v G L«(fi,M'^'), -> v strongly m L'^{n x Y,R''), 

(b) Let be a bounded sequence in Li{n,M.^) such that ^ V strongly in 
L'J{n,R^), then 

Tri{vn)^v, strongly in L'^{Q X Y,R''). 

(c) For every relatively weakly compact sequence Vj^ in L'^{n,M.'^), the corre- 
sponding Tjjivri) is relatively weakly compact in L^(ri x y, R'^). Further- 
more, if 

%j{vr,) in L«(rj X r,M'^'), 

then 



1 



vdy in L'^{n,R^). 



Proof See [H Proposition 2.9]. □ 

We note that the strong/weak convergence of Tri{vrj) does not imply, unless 
has the measure for every 77, the strong/weak convergence of w,,, since the 
information concerning the behavior of f,, on A,j is missing. 

Next results present some properties of the restriction of the unfolding op- 
erator to the space W^^'^{n,R''). 
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Proposition 3.4. Let q belong to ]1, oo[. 

(a) Suppose that Vjj G W^''^{ft,M.'^) is bounded in L^(J7,]R'^) and satisfies 

v\\^Vn\U<C. 

Then, there exists a subsequence and v e LP{n,Wp^1{Y,M.^)) such that 

(b) Let Vjj converge weakly in W^''^{fl,R'') to v. Then 

Proof. See [HI Corollary 3.2, Corollary 3.3]. □ 

Proposition 3.5. Let q belong to ]1 , oo[. Let converge weakly in W^''^{il,W^) 
to some V. Then, up to a subsequence, there exists some v G L'^{Vl, Wpi^iY, R*^)) 
such that 

TriiVv^) ^Vv + Vyi) in L'i{VtxY,M}'). 

Proof. See [131 Theorem 3.5, (i)]. □ 

The last proposition can be generalized to ]R'^)-spaces with m > 1. 

Proposition 3.6. Let q belong to ]1,cxd[ and m > 1. Let w,, converge weakly 
in W'^''^(yi,M.^) to some v. Then, up to a subsequence, there exists some v G 
L«(17,M^p™^«(y,R*=)) such that 

TriiD^Uri) D^v in VF'"-''''(F,R'=)) for \l\ < m - 1, 

TrjiD^v,^) D^v + D'yV in L^iQ x Y,m'') for \l\ = m 

Proof See [l3l Theorem 3.6]. □ 

We note that the periodic unfolding method described above is an alternative 
to the two-scale convergence method introduced in [46 and further developed 
in [B]. More precisely, the two-scale convergence of a bounded sequence Vrj in 
LP{^,M.'') is equivalent to the weak convergence of the corresponding unfolded 
sequence Tr,{v^) in LP{n x Y^R^) (see [TF, Proposition 2.14]). 

For a multi- valued function h G M(r2, M*^, a, mjl we define the unfolding 
operator as follows. 

Definition 3.7. Let Y be a periodicity cell, rj be a positive number and a map 
h G M(f],R'=,p,a,m). The unfolding operator %{h) : n x Y x ^ 2^" is 
defined by 

jh(7]\^] +7]y,v), a.e. (x, y) G a, X y, G R^ 
\\v\P-''^v, a.e. G X r. w G M'^. 



'The class of functions h G M(r2, K**, a, m) is defined in Definition l5.ll 
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Obviously, by its definition the unfolding operator of a multi- valued function 
from M(f],R'=,a,m) belongs to the set M(f7 x Y,R'',a,m). 

We note that the periodic unfolding method described above is an alternative 
to the two-scale convergence method introduced in [33] and further developed 
in [6]. More precisely, the two-scale convergence of a bounded sequence Vri in 
Lp{Q,'EJ') is equivalent to the weak convergence of the corresponding unfolded 
sequence Tr,iv^) in LP{n X y,]R'=) (see Proposition 2.14]). 

Homogenization of the linear elasticity problem. In this section we ap- 
ply the periodic unfolding method to the homogenization of linear elasticity 
system^l^ with periodically highly oscillating coefficients (see [HJ for properties 
of periodically oscillating functions). We show the strong convergence of the 
unfolded sequence of the gradients of the solutions of linear elasticity problem 
(see Theorem 13.81 below). The proof of the mentioned result applied to an ellip- 
tic partial differential equation is performed in |13| and can be carried over to 
linear elasticity systems without significant modifications. Therefore, we sketch 
here only the proof in [13] adopted to our needs. 

In linear elasticity theory it is well known (see [55] Theorem 4.2]) that a 
Dirichlet boundary value problem formed by the equations 



to given b G iJ^^(57,M^) and G L^{il,S^) has a unique weak solution 
(w^,tT^) e i/o(^'-'^'^) ^ L^{i^,S^). We require that converges to Eq strongly 
in L-^{fl,S^) as ?7 ^> 0. The following result holds. 



Theorem 3.8. There exist uq G H^{n,M.^), ctq G L'^{ft x Y,S^) and ui G 
(y,M3)) such that 



an(x) 
Unix) 



b{x), X fl, 

C[x/ry](sym(Va;M,,(a;)) -e^(x)), 

0, X G dn, 



X eO, 



(15) 
(16) 
(17) 



It,, ^ Mo mi7o(r2,R^), 
%{un) uo m L^{n,H^^,{Y,R^)), 
%{Vur,) Vuq + VyUi in L'^{n x Y,W^ 
7;,(cr„) ^(70 in L^in X Y,S^), 



) 



(18) 



(19) 



(20) 
(21) 



and (uo,(7o,ui) is the unique solution of the homogenized system: 



divy CTo(x, y) 
CTo(a;,y) 




y 




(25) 



uo(x) 



= 0, 



X G on. 



(26) 



survey on other applications of the method can be found in |13| . 
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Moreover, the following convergences hold 



lim / C [-/r;] sym(V2;U„) syni(Va;M„)(ia; 



= / C[y] sym (\7 uq + VyUij sym (\/uo + VyUi^dxdy, (27) 
lim / I syni(Va;M„)pda; = 0, (28) 



%{Vu^) ^Vuo + VyUi in L^VLxY.W'), (29) 
7;,(ct,,) ^ CTo «n -/-^(f^ X y,^^). (30) 
Proof. See H5| Theorem 4.1]. □ 

4 Unfolding the Curl Curl-operator 

Our method is based on the Helmhohz-Weyl decomposition for vector fields in 
general L'^-spaces over a domain £7. It turns out (see [301 Theorem 2.1(2)]) that 
the following theorem holds. 

Theorem 4.1. Let 1 < q< oo. Every v G L'^(ri,M"^) can he uniquely decompose 
as 

V ^ h + G\\r\w + Vz, (31) 

where h e Xl^^{n,M?), w e V;i{n,M?) and z e W^^'i{n,M?), and the triple 
{h, w, z) satisfies the inequality 

\\h\U + \\w\\u,si + Mi.g.n < C\\v\\„ (32) 

where C is a constant depending on fl and q. If there is another triple of 
functions {h,w,z) such that v can be written in the form 

V ^ h + Curl w + Vz, 

with h e Xl^^{n,R^), w e V^{n,R^) and ~z e W^''i{^,W'), then it holds 

h — h, Curlw — CurlTi, Vz = Vz. 

Remark 4.2. If L denotes the dimension of Vhari^i^^), i-e. dlTaVhar{^^R^) = 
L, and {(/>!, ...,02} is abasisof y,,ar(f^,K^), then it holds F«(r2,M3) c W^^''{n,R^) 
with the estimate 

L 

\\v\U + l|Vt;||, < C{\\ Cm\v\\, + II div^ll, + J2 l("i<^OI) 

1=1 

for all G F«(f7,]R3)^ ^here C = C{n,q) ([501 Theorem 2.4(2)]). The proof 
of the inequality with X^iLi l('^:<^i)l replaced by ||w||g is performed in [301 
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Lemma 4.5] (for q = 2 it can be found in [THl Theorem VII. 6.1]). If we assume 
that the boundary dQ has L + 1 smooth connected components ro,ri, ...jF^ 
such that Fi, he inside Fo with F,; n Fj = for i ^ j and 

then it holds ([30l Appendix A]) 

If the function v in (j3ip is more regular, then the function w can be chosen 
from a better space as the next theorem shows. 

Theorem 4.3. Let 1 < q < oo. Assume that decomposition i31\) holds. If, 
additionally v e Z^^^^{il,R^), then w in if^TI) can be chosen from W^''^{V,,R^)n 
V^(r2,]R^) satisfying the estimate 

Ikll3,g,n <C(||Curli;||i,g,o + ||«||g), (33) 
where C is a constant depending on D, and q. 

Proof For v G L^^^^j(17, R^) this resuU is proved in For v e Zl^^.^{n,M.^) 
the proof runs the same lines. We repeat them. 

As it is shown in [301 Lemma 4.2(2)], we can choose the function w e 
V^{il,M.^) satisfying the equation 

(Curlw, Curl V')ji = (w. Curl ?/.)o, ior all G Vf {n,R^) (34) 

with the estimate 



\Mi,,,n < C\\v\\g, (35) 

where C depends only on and q. Since divw = in f2 and v £ Z^^^^{i},R^), 
it follows from that —Aw = Curlw in the sense of distributions, and we 
may regard w as a weak solution of the following boundary value problem 

— Aw — Curlu, in f2, (36) 
divti; = 0, on dfl, (37) 

w-n = 0, on dn. (38) 

Since Curlw e W^''^{il,R^), it follows from [SUl Lemma 4.3(1)] and the classical 
theory of Agmon, Douglas and Nirenberg W that the solution w of the homo- 
geneous boundary value problem (l36)) belongs to W^''^{i},R^) and the estimate 

IHU.g.n < Ci\\Cnvlv\\i,,^n + \\w\\g), (39) 

is valid with the constant C dependent of ft and q. Due to the estimate 
implies ([55)1 . This completes the proof. □ 

Now we can state the main result of this section. 
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Theorem 4.4. Let \ < q < oo. Suppose that sequence is weakly compact in 
■^Curi(^' ^'^)- Th^''^ there exist 

G Z«^^i(f2,R3), Wo G i^C^^ X F,R3) and 



such that 

Vri^v in Z^c^jn,R^), (40) 

%{v^) ^ «o m L'^i^ X y,R3), (41) 

r^(Curlv^) Curlw in L'J{n,W^f^{Y,R^)), (42) 



7;,(CurlCurlu^) Curl Curl + Curlj, Curl^, ui , in L^in x y,R3).(43) 
Moreover, v{x) — Jy voix,y)dy. 

Proof. Convergence (|4ip and the last statement of the theorem follow from 
Proposition I3.3f c) . Convergence is obvious. Next, we prove convergences 
dH]) and (031). According to Theorem Ol there exist /i^ G Xl^^{n,R^), Wr, G 
V;«(fi,M3) and z,, G W^i'«(rj,R3) satisfying the inequality 

WKWq + \Wr,\\l^q.n + Ikr,||l,g,fl < (44) 

with the constant C independent of 77, and such that 

Vri = + Curlw^ + Vz^. (45) 
Moreover, due to Theorem 14. 3[ in enjoys the inequality 

< C{\\ Curl«^||i,,,n + ||«^||,) (46) 

with the constant C independent of rj. Therefore, the weak compactness of in 
Z^^^^{n,R^) and ge]) imply that Wr, is weakly compact in W^''i{n,R^). Thus, 
in virtue of Proposition 13.61 we conclude that there exist 

w (£W^''^{n,R^) and wi e L\n,W^^^^{Y,R^)) 

such that 

%,{D^w„) ^ D^w in L'i{n,W^-''''{Y,R'')) for \l\ < 2, 
Tr^iD^w,,) ^ D^w + D^wi in L'^{Q. x Y,R^) for \l\ = 3. 

Since Curlw,, = Curlu;,, and Curlw = Curlw, we get that 

r^(Curl?;^) ^ Curlv in L''{n,Wl;!^{Y,R'^)), 

7;,(CurlCurlw,,) Curl Curl w + Curl^, Curl^, wi , in L'^{n x Y,R^). 
The proof of Theorem 23] is complete. □ 
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5 Homogenized system of equations 

Main result. First, we define a class of maximal monotone functions we deal 
with in this work. 

Definition 5.1. Form € Li(n,E), a e IR+ and q > 1, M{n,R'',q,a,m) is the 
set of multi-valued functions h : V, x MJ^ 2^ with the following properties 

• V i-^ h(x,v) is maximal monotone for almost all x ^ fl, 

• the mapping x i~> j\{x,v) : 17 — )• M*^ is measurable for all A > 0, where 
j\{x,v) is the inverse of v t-^ v + Xh{x,v), 

• for a.e. x G O and every v* € h{x,v) 

ai^ + ^^j<{v,v*)+m{x), (47) 

where l/q + 1 /q* ~ 1 . 

Remark 5.2. We note that the condition (|T7|) is equivalent to the following two 
inequalities 

\v*f < mi{x) + ai\v\'^ , (48) 
(u, V*) > m2{x) + a2\v\'^, (49) 

for a.e. x £ fl and every v* £ h{x,v) and with suitable functions mi,ni2 G 
L-'^(f2,R) and numbers ai,a2 G R+- 

Remark 5.3. Visco-plasticity is typically included in the former conditions by 
choosing the function g to be in Norton-HofF form, i.e. 

g(I]) = [|I]|-ay];^, i:eM\ 

where Cy is the flow stress and r is some parameter together with := 
max(x, 0). If (? : i-> then the flow is called irrotational (no plastic spin). 

The main properties of the class M(il, M'^, g, a, m) are collected in the fol- 
lowing proposition. 

Proposition 5.4. Let % be a canonical extension of a function /i : M'^ — > 2^ , 
which belongs to M(r2, R*^, q, a, m). Then Ti is maximal monotone, surjective 
and D{n) = LP{n,R^). 

Proof See CoroUary 2.15 in [TC]. □ 

Next, we define the notion of strong solutions for the initial boundary value 
problem © - ([TO)) . 

Definition 5.5. (Strong solutions) A function {uri,(Jri,Pri) such that 

is called a strong solution of the initial boundary value problem (0) - ilO\) . if for 
every t £ [0, Tg] the function (urjit), ari{t)) is a weak solution of the boundary 
value problem U5\) - Jj?] ) with ip = symp^(i) and b — b{t), the evolution 
inclusion ^ and the initial condition ^ are satisfied pointwise. 
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Next, we state the existence result (see '45' ) . 

Theorem 5.6. Suppose that l<q*<2<q< oo. Assume that fl C M.'^ is a 
sliceable domain with a -boundary, C G L°°{n,S^) and Ci G L°°(n,M). Let 
the functions b G W^i'«(0, Te; L9(f7, M^)) fee given and g G M{n,M^,q,a,m). 
Suppose that for a.e. x E D, the relations 

p'-°\x)=0 and G g(a;/?7,CT(°)(a;)) (50) 

hold, where the function a^^^ G L^{^l,S^) is determined by equations m5\) - JiTp 
for ip = and b — 6(0). Then there exists a solution (m,,, cr^,p,,) of the initial 
boundary value problem (0) - ilO\) . 

Now we can formulate the main result of this work. 

Theorem 5.7. Suppose that all assumptions of Theorem \5.6\ are fulfilled. Then 
there exists 

uoeH'{0,T,;H^{n,R'')), uieH\0,T,;L^{n,Hl^,{Y,R^))), 
ao G L°^{0,T,;L\n x Y,S^)), pa G H\0,T,;L^in x Y,M^)), 
p G H\0, Te; L^n, n 2.2(0, Te; Z^^uA^, 

pi G L\0,T,;L^in,W^'f iY,M^))) with divyPi=0, 

and 

such that 

urj^uo in H\0,T,;H^in,M.^)), (51) 
p^^p in H\0,T,;L^{n,M''))nL\0,T,;Zl^^,{n,M^)), (52) 
7;,(Vu^) ^ Vwo + VyMi in H\0,T^;L^{n xY,W^)), (53) 
a„^ao ^n i°°(0,r,;L2(r!,53)), (54) 
%{a„)^ao m L°^{0,T,;L\nxY,S^)), (55) 
7;,(Pr,)^Po m L\0,T,;L^nxY,M^)), (56) 
7;,(atp^) ^atpo in L2(o,re;L2(l] x r,^^))^ (57^ 

and 

r^(Curlp,) Curlp ^n ^^(o, T^; ^^(l], ■^^))), (58) 

7^(devsymp,,) ^ devsympo in L^{D,t^ xY,M^), (59) 
7;,(CurlCurlp,,) m ^^(fi^^ x F, X^)^ (gQ) 

r,(i]!;") - 4'" *n L''(f]T. xy,x3)^ (gi) 

w/iere 

p := Curl Curlp + Curlj, Curlj^ pi, 
Sq" := o-Q - Ci[y] devsympo - C2P, 
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and (uq, Ui, cr, aQ,p,pQ,pi) is a solution of the following system of equations: 



~dw^ a{x,t) = b{x,t), (62) 
~divyao{x,y,t) = 0, (63) 

cro{x,y,t) = C[y]{sym(VxUQ(x,t) +WyUi{x,y,t) - po{x,y,t)))l64:) 

dtPoix,y,t) e 5(2;,4'"(a;,y,t)), (65) 

which holds for (x, y,t)^nxM.^x [0, Tf,], and the initial condition and boundary 
condition 

p{x,0) = p'-°Hx), xen, (66) 
p{x,t)xn{x) = 0, {x,t) Ednx[0,Te), (67) 

uo(a;,i) = 0, {x,t) ednx[0,Te)., (68) 



T/ie functions a and p are related to ctq and po in the following ways 

a{x, y cro(a;, y, t)dy, p{x, = J Po{x, y, t)dy. 

The proof of Theorem 15. 71 is divided into two parts. In the next lemma we 
derive the uniform estimates for (u^, (7^,p^) and then, based on these estimates, 
we show the convergence resuh. 

5.1 Uniform estimates 

First, we show that the sequence of solutions (u,,,(7^,p^) is weakly compact. 

Lemma 5.8. Let all assumptions of Theorem \5.7\ be satisfied. Then the sequence 
of solutions (m^,(t^) is weakly compact in _ff^(0, Tg; i/g (fJ, R'') x L'^{fl,S^)) and 
Pn IS weakly compact in H^{0,Te; L^ifl, M^)) D L^{0,Te, Z^^JQ, M^)). 

Proof. To prove the lemma we recall the basic steps in the proof of the existence 
result (Theorem 15. 61) . For more details the reader is referred to 35] ■ The time- 
discretized problem for ([5]) - ([TU)) is introduced as follows: 
Let us fix any to e N and set 

rp -1 />7lh 

/i:=T5^, <™:=0&'^:=-r / b{s)ds e L'^in,R^), n = l,...,2™. 

^ "J(n-l)/i 

Then we are looking for functions G H^{Q,'R^), ct^^ e L'^{fl,S^) and 

PjJ.rn Zl^^^{n,M^) with p!^^^{x) e s[(3) for a.e. x e f2 and 

<™ - C^ [x/77] dev symp;;_„ - -pl^^ - C2 Curl Curlp^J^^ e (17, M^) 



solving the following problem 

-div,a;;,„(x) = b^,{x), (69) 

al^{x) = C[x/rf^{syu,{V^ul„^{x)-pl„^{x))) (70) 

<^(^)-pr^^(^) ^ ,(,/,, E;j"„(x)), (71) 
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together with the boundary conditions 

pl^{x)xn{x) = 0, xedn, (72) 

ul^ix) = 0, a; ear!. (73) 

Such functions (uJJ^m, crJ^_m,p5J_„i) exist and satisfy the following estimate 



^{m'^'<^n.J\l + «i|l devsymp^^„||2 + ^Wp'^^Jj + C2II CmW^^J\l 

„n-l 



+hCY: \K:m\\l + f If , < + / mix)dx (74) 

I 

+hcY,(\rj\l + \W^-bl;')/h\[t 



for any fixed / G [1, 2™], where (here 



2C(") := |lBi/2a(°)|| 



and C, C are some positive constants independent of 77 (see [J^ for details) . To 
proceed further we introduce the Rothe approximation functions. 

Rothe approximation functions: For any family {Cm}n=o,...,2m of functions 
in a reflexive Banach space X, we define the piecewise ajfine interpolant e 
C([0,Te],X) by 

U{t) ■■= {{-{n- I)] Cn + - for {n~l)h<t< nh (75) 



and the piecewise constant interpolant g L°°(0,Te;X) by 

U{t) C for {n-l)h<t< nh, n^l, 2'", and e°„. (76) 

For the further analysis we recall the following property of and S^m'- 

, (77) 

where is formally extended to t < by and 1 < q < 00 (see |51jV 
Now, from (I74p we get immediately that 



(Olli2 + "i||dev Sym Prj.m ml + -\\Pv.rn{t)\\l + C2II Curlp^,„(t)||2 

+2(7 (||9tP,,™||^:,o,(o,T.) + PJjril^,^,x(o,T.)) (78) 
< 2C(°) + 2||m||i,o + 2C||&f^,„(„^^^^^,(„ .,3)), 

where C is some other constant independent of ry. In [45 it is shown that the 
Rothe approximation functions {uri^m,cfn,m,Pri,m) and {ur^^m^o'ri^rmPri.rn) con- 
verge to the same limit {u^i, ari,Pri)- Due to the lower semi-continuity of the 
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norm and (|78l) this convergence is uniform with respect to 77. Therefore, esti- 
mate (175]) provides that 

{an}ri is uniformly bounded in L°°(0, T^; £^(1], 5^)), (79) 
{devsymp j^ is uniformly bounded in L°° {0,Te; L^{n, M^)), (80) 

{Curlp^}^ is uniformly bounded in L°"{0,T^; L^{n,M^)), (81) 

{p,,},, is uniformly bounded in W^''^' (0, Te; L-?* (17, TW^))^ (32) 

{Sjj"},, is uniformly bounded in (17t, , M^). (83) 
Furthermore, from estimates ([75]) - we obtain easily that 

{u^}^ is uniformly bounded in L^(0, T^; H^{n, R^)), (84) 

{p^},, is uniformly bounded in L^{0, T^, ^Suri(^: -^^))- (85) 

Additional regularity of discrete solutions. In order to get the additional 
a'priori estimates, we extend the function 6 to i < by setting b{t) — 6(0). The 
extended function b is in the space W^'P{-2h,Te; W~^'P{n,R^)). Then, we set 
= b^i^ := b{0). Let us further set 

where ^/^ : L^(f2,A^'^) — > 2^'''^^''''^'^)) denotes the canonical extensions of g(x/r/, •) : 

2='(3) The assumption §U\) implies that p^^^ = 0. Next, we define func- 
tions (u^^m,cr~^) and (M°_m)^r;,m) solutions of the linear elasticity problem 
(pB - (dZl) to the data & = 7 = 0, = and & = 6^, 7 = 0, Ep = 0, 
respectively. Obviously, the following estimate holds 





7/0 — 7, 1 




?7,m ?7,m 






h 


2 




J 



where C is some positive constant independent of m and 77. Taking now the 
incremental ratio of ((7T|) for n = 1, 2™, we obtaird 

Let us now muhiply the last identity by - ,„)//i. Then using the 

monotonicity of Qrj we obtain that 

^(rtp;;,™ -rtp;-\rtp;; „) + (rtp;;^„ -rtp;-\Ci devsym(rtp;_„))^ 

+ (rtp;^„ -rtp,^;„\ 6*2 Curl Curl(rtp;;^„J)^ < ( rtp;;^„, - rtp,^-i, rt C7;%J^. 
With (|69|) and (j70l) the previus inequality can be rewritten as follows 

:;^(rtp;;,„ -rtp;;;;,\rtp;;,„) + (rtp;;,„ -rtp;\„\Ci devsym(rtp;;„,))^ 

*For sake of simplicity we use the following notation rt</>5^ := {4>m ~ 4>m^^)/h, where 
'/'m' 'Ami is any family of functions. 



23 



+ {i-tpl^-i-tp--^,C2GuT\Gm\{rtplJ)^^ + ( rt - rt C'^-t 

As in the proof of (j74p . multiplying the last inequality by h and summing with 
respect to n from 1 to I for any fixed / G [1, 2™] we get the estimate 

^11 rtp!,,„|l^ + toill devsymrtp(,_„|l2 + h\\M'^' vt a^^^jH + hC^W Curlrtp^^^H^ 
I 

< 2h&'^ + 2hY^ (rt<„ - rt<;„\rt6;5,)^, (87) 
where now C(°) denotes 

2CW :^||BV2rt<„||. 

We note that (j86p yields the uniform boundness of C^^-* with respect to m. Now, 
using Young's inequality with e > in (j87p and then summing the resulting 
inequality for I = 1, 2™ we derive the inequality 

^ ll^tPmllz.nT, +"1 l|devsym(atp,„)||2_o^^ + C2 ||Curl (atp™)||2 (88) 
+C||at(T„||^^^^^ < C,\\dtb^\\i^^^+2e\\dtu^\\in^^, 

where is some positive constant independent of m and 77. Using now inequal- 
ity (HI), the condition dtPm{x,t) € sl(3) for a.e. (x,t) S f^Te, and the ellipticity 
theory of linear systems we obtain that 

^ \\dtPrn\\ln^^ +C,in) \\dtPm\\ln^^+C\\dta„Ml^^^ < C^dtbmWln^^, (89) 

where Ce(i7) is some further positive constant independent of m and ry. Since 
bm is uniformly bounded in W^''^{nT^,S^), estimates (I55|) and ([55]) imply 

{devsymcJtp^}^ is uniformly bounded in L^(0, Tg; L^(r2, A^^)), (90) 

{9(0-,,}^ is uniformly bounded in L^{0, T^; L^{n, M^)), (91) 

{Curiatp^}^ is uniformly bounded in L^{0,Te; L^{n,M^)), (92) 

{p^},, is uniformly bounded in H^{0,Te; L^^.^in, M^)). (93) 

The proof of the lemma is complete. □ 

5.2 Proof of Theorem EH] 

Now, we can prove Theorem 15. 71 

Proof. Due to Lemma [5.81 we have that the sequence of solutions {uri,<Tri) is 
weakly compact in i?^(0, Tg; iJp (il, M"^) x L^{n,S^)) and the sequence p,j is 
weakly compact in H^{0,Te; L^{n, M^)) n L^{0,Te; Z^^^^in, M^)). Thus, by 
Proposition 13. 3[ Proposition 13.51 and Theorem 14. 4[ the uniform estimates ([7^ 
- (|93p yield that there exist functions uq, ui, a, uo, p, po and pi with the 
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prescribed regularities in Theorem 15 .61 such that the convergences in ([ST|) - (pT|) 
hold. Note that - 1^ give the equation ([M]), i.e 

cro{x,y,t) = C[y]{sym{WxU(){x,t) + WyUi{x,y,t) ~ po{x,y,t))), a.e. (94) 

By Proposition 13.31 the weak-star limit a of ct^ in L°°{0,Te; L^{il,S^)) and 
the weak limit p of p,, in L^{0,Te; L^{il, A4^)) are related to (Tq and po in the 
following ways 

a{x,t) = J ao{x,y,t)dy, p{x,t) = J po{x,y,t)dy. 

Now, as in 171, we consider any e C|^(r2, M"^). Then, by the weak convergence 
of (7^, the passage to the weak limit in ^ yields 



{a{x,t),V(t){x))dx ^ / {h{x,t),(t){x))dx, (95) 

i.e div^ CT = 6 in the sense of distributions. Next, define 4>n{x) — ri(j){x)'ip{x/r]), 
where (f) £ C^{n,M.^) and ip € C^^{Y,R^). Then, one obtains that 

Therefore, since 0^ has a compact support, 

{%j{<T^{t)),%{\7(l)^))dxdy^ [ {h{t),(j)^)dx. (96) 



inxY Jn 
The passage to the limit in ([M|) leads to 

{(Ta{x, y, t), (j){x)Vyip{y))dxdy = 0. 

nxY 

Thus, in virtue of the arbitrariness of (j), one can conclude that 

iaoix,y,t),Vyij{y))dxdy^O. (97) 

nxY 

i.e divy (Tq^x, - ,1) — in the sense of distributions. 

Next, let %{g,^) : LP{VL x r,M^) ^ 2^''(^^><^'K") and ; LP(1],R^) ^ 
2L''(n,E") (igj^o^g the canonical extensions of %{g^){x,y) : 2^" and 

5(j/) : — >■ 2"^ , respectively. Here, g{y) is the pointwise limit graph of the 
convergent sequence of graphs Tr]{gri){x, y). The existence of the limit graph for 

%f{grj){x, y) guaranteed by Theorem 12. 71 Indeed, the resolvent converges 
pointwise to the resolvent , what follows from the periodicity of the mapping 
y ~^ giUj 2) : y — )■ 2 and the simple computations: 

jI^^'^\x, y, z) = TM'){x, ^) = Jiiv^ 
for a.e. (x, y) G x F and every z G M^. Thus, by Theorem 12.71 we get that 

rM{^^y)^9{y) (98) 
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holds for a.e. {x,y) € fl x Y. Since g,, 6 7W(n,IR^,p,a,m), by Defini- 
tion 13.71 of the unfolding operator for a multi-valued function it follows that 
%j{g,f) & M{fl X y, E^,p, a, m). Therefore, due to this and convergence (^S)). 
by Propositon l5.4f b) we obtain that 

%ig^)>^g. (99) 



To prove that the limit functions ((To,po) satisfy (|65p. we apply Theorem 12.61 
Since the graph convergence is already established, we show that condition ([1 
is fulfilled. Using equations ^ and we successfully compute that 



1 

W 
1 

W\ 



nxY 

{Tr,{dt{e{Vur,{t)) - C-^a^n{t))) ,Tn{ar,{t)))dxdy 



lin (4.\ I 



= / {h{t),dtUr,{t))dx- / {b{t),dtUr,{t))dx 

Jn J A,, 

-j^ I {%{dtC-'a,m,rMt)))dxdy 
inxY 



1 

m 
1 

W\ 



nxY 



(CiTriidt dev sym Pr,{t)),%j (dev sym Pr,{t))) dxdy 



{C2%, (dt Curl p^ (t)), Tr, (Curl p,, {t)))dxdy. 
nxY 



Integrating the last identity over (0,i) and using the integration- by-parts for- 
mula we get that 



(r^ (5tp^ (t ) ) , (Sl)" (t ) ) ) o X y (100) 
{b{t),dtu^{t))ndt- / {h{t),dtu^{t))^Jt 



1 



-\mB''^a,m\lnxY + \\\Tn{B"'<y,m\\lnxY 
-^ry'r,(devsymp,(0)||loxy - ^r2'/'7;,(Curlp,(i))||2 

where ;B = C^^. Due to the uniform boundncss of {9tu^} in L^(0, Tg; i7Q(i7, 
we easily obtain that 



lim / {h{t),dtUr,{t))t,Jt^Q. 

Moreover, since (T,,(0) solves the linear elasticity problem ([T5|) - ([T7]) with e,, = 
and b = b{t), by Theorem l3.8l we can conclude that 7^(S^/^cr,,(0)) converges to 
B^^^ao{0) strongly in L^{D. x Y,S^). Thus, by the lower semi-continuity of the 
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norm the passing to the limit in (|100p yields 
1 /■* 

limsup-- / {%{dtP^it)),%i^'^''it)))n^Ydt 

n— >oo I I Jo 

-i||Ci devsympoWll2,nxy " ^l|C2 CmlpoWllloxy, 



or 



1 



Y 



limsup — / {%{dtP^{t)),%{E\]''{t)))nxYdt 



n—>-oo 

< ^ {b{t), dtUoit))ndt - ^ / (atC-Vo(t), ao(t))j^^ydt (101) 
t 



1 /■ 

- 1^ y (9t dev sym po ) , Ci dev sym po (i ) ) 



We note that ((Ml) and ^ imply 

[ {b{t),dtUo{t))dx^-^ [ {ao{t),dts{^uo{t) + \/yUi{t)))dxdy. (102) 
And, since for almost all {x,y,t) E il x Y x (0, T^) one has 

{dt dev sympo(a;, y, i), Ci [y] devsympo(a;, y, t)) = (dtPoix, y, t), Ci [y] devsympo(a:^ 
and that for almost all t G (0,Te) 

(9tCurlpo(t),C2Curlpo(i))f,^^ = (9tPo(i), C2 CurlCurlpo(0)axy' 
the relations (ITOTI) and (fT02l) together with dMl) yield 

limsup— / {%{dtP^{t)),%{Y}:^mnxYdt 

n— >-oo l-f I Jo 



/O 



{dtPo{t),T}^{t))^^^dt. (103) 



In virtue of convergence (|99l) and inequality ([TOS)) . Theorem [2Jl vields that 
[Sj;"(x,y,i),atPo(a;,y,0] eGr5(y) 

or, equivalently, that 

dtPo{x,y,t) e g{y,Y.a'^{x,y,t)). 

The initial and boundary conditions (j66l) - (j68p for the limit functions uq and 
p are easily obtained from the weak compactness of and Prj in the spaces 
H^O,Te;H^in,R^)) and iJ^O, L2(f7, TW^)) n /.^(o, T^; Z2^,i(f}, X^))^ 

spectively. Therefore, summarizing everything done above, we conclude that 
the functions (wq, mi, cr, uoiPiPo 7 Pi) satisfy the homogenized initial-boundary 
value problem formed by the equations/inequalities ([5^ - □ 
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